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Abstract. In this article we study the well-posedness (uniqueness and exis- 
tence of solutions) of nonlinear elliptic Partial Differential Equations (PDEs) 
on a finite graph. These results are obtained using the discrete comparison 
principle and connectivity properties of the graph. This work is in the spirit 
of the theory of viscosity solutions for partial differential equations. 

The equations include the graph Laplacian, the p-Laplacian, the Infinity 
Laplacian, and the Eikonal operator on the graph. 



1. Introduction 

In this article we consider existence and uniqueness of solutions to nonlinear 
elliptic partial differential equations on a finite graph. The uniqueness results are 
based on appropriate versions of the the discrete comparison principle, which for 
some equations depends on fine connectivity properties of the graph. 

1.1. Context and motivation. The graph Laplacian has been of interest since 
Birkhoff studied it in the nineteenth century. Finite difference discretizations of 
nonlinear elliptic Partial Differential Equations (PDEs) have been used extensively 
in Image Processing [Sap06] . Modern applications of PDEs on graphs include ma- 
chine learning, clustering, and social networks [GP10]. 

By their very definition, PDEs on graphs use only local information: locality is a 
requirement for problems on distributed networks [OT11]. Very large graphs chal- 
lenge the primarily combinatorial tools originally designed to study them [Lov09]. 
One popular measure of tractability is that algorithms run in polynomial time. For 
our purposes, that is not sufficient. The PDEs of the types we study here can be 
solved quickly (in log-linear time) on grids [Tsi95, TCOZ03, FOll]. This kind of 
fast solvers may also be available for this class of equations on general graphs. 

Nonlinear PDEs on graphs arise naturally in stochastic control theory, when the 
state space is discrete [Ber05]. Numerical methods for stochastic control problems 
can be shown to converge using probabilistic methods [KD01] or using the viscosity 
solutions theory [BS91]. In general, wide stencil schemes are needed to represent 
linear equations by positive difference schemes [MW53]. Most discretizations as- 
sume a compact stencil, which leads to specific requirements on the structure of 
the equation, such as diagonal dominance [BS91]. Extending these schemes to gen- 
eral problems requires the use of wide stencils [BZ03]. More recently, Lions and 
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Lasry have also extended the H JB theory to Mean Field Games [LL06] . Mean Field 
Games have been posed on graphs [Guell,GMS]. 

Stochastic control theory in the continuous setting leads to Hamilton- Jacobi- 
Bellman (HJB) equations [FS06]. These equations motivated the theory of viscos- 
ity solutions. The theory of nonlinear (and possibly degenerate) elliptic PDEs in 
continuous space is now well understood [CIL92] . In general, discretization of these 
PDEs using the finite difference method may not converge. However, if the finite 
difference schemes also obey a comparison principle, then a convergence proof is 
available [BS91]. These finite difference schemes, which are called elliptic [Obe06] 
can be characterized by a structure condition. Elliptic schemes generalize upwind 
schemes for first order equations and schemes of positive type for second order 
equations [MW53]. The PDEs studied here on general graphs coincide with ellip- 
tic finite difference schemes when the graph is a wide stencil finite difference grid. 
Finite differences are discussed in Section 7. 

The Infinity Laplacian PDE [ACJ04] has been well studied in both the discrete 
and continous settings. A convergent discretization of the Infinity Laplace equations 
was presented in [Obe05] (see also [LGA98]). This discretization is reinterpreted 
here as a PDE on a graph. This equation has also been interpreted using random 
tug-of-war games, in both the continuous and the discrete setting [PSSW09]. Vari- 



ants of the Infinity Laplace equation can be posed on graph (see (3.16) below), these 
also arise as finite difference schemes for p-harmonic functions [Obel2]. The unique- 
ness for Infinity Laplace equation follows from the uniqueness of the finite difference 
scheme, and from the comparison with cones property [AS10, AS12, APSS1109]. 
These results have been extended to more general equations in [ACJS11] and to 
equations with drift in [ASS11]. The existence and uniqueness of the solution of p- 
Laplacian and the connection with the game interpretation has been futher studied 
in [PS08] and [MPR10]. 

An elliptic discretization for the equation for motion of level sets by mean cur- 
vature was presented in [Obe04] . An earlier work [CDK95] presented a convergent 
method, but this was semi-discrete: although the method was implemented on a 
grid, the scheme was presented in the continuous case. The scheme in [Obe04] 
involved the median (see (3.5) below): the use of wide stencils was required to fully 
discretize the equation onto a grid. The game interpretation of motion by mean 
curvature was presented in [KS10]. The game involved a formula similar to the one 
in [CDK95]. The median scheme and the game scheme are both consistent, so they 
yield very similar results in the continuous setting. This game formulation does 
not have a natural graph interpretation because it requires a notion of direction 1 , 
which is not available on a graph. The median scheme is defined on a graph below. 
However, there is no reason to expect in general that solutions of the resulting flow 
on the graph respect the same properties (such as the shrinkage of level set curves) 
as do solutions of the continuous PDEs. 

1.2. Results of the article. It is natural to ask whether there are general condi- 
tions on the PDE which lead to existence and uniqueness results. In this article, 
we establish well-posedness (uniqueness and existence of solutions) results for non- 
linear elliptic PDEs using structure conditions on the equations and connectivity 
properties of the graph. We do not rely on linearity, variational structures, game 

^One player choses a direction vector, the other player choose whether to move in the direction 
of the vector or the opposite direction. 
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interpretations, or optimal control interpretations. Instead, we use the comparison 
principle as the basic tool. 

We begin with the two simplest structure conditions. By analogy with [CIL92] 
these correspond to proper and uniformly elliptic equations. Existence and unique- 
ness results for proper equations were established in [Obc06]. Our first result here 
is to define uniformly elliptic equations and to prove well-posedness. In the uni- 
formly elliptic case, the uniqueness result uses the idea of marching to the boundary, 
which can be found in the early paper of Motzkin-Wasow on linear elliptic finite 
difference schemes [MW53]. This proof readily generalizes to the nonlinear case. 
Well-posedness can fail for elliptic PDEs on graphs, as examples below show. Our 
results are specific enough to avoid these examples, while still being general enough 
to consider a wide class of operators. 

The next class of operators we consider are generalizations of the Eikonal and 
Infinity Laplace equations. These are degenerate (neither proper nor uniformly 
elliptic) equations which depend on the maximum and the minimum of the neigh- 
boring values. The uniqueness proof for the generalized Infinity Laplace equations 
is a variation of the proof in [LGA98], see also [AS10] and [AS 12]. Another proof 
based on martingales is in [SvilO]. 

Remark 1. We do not discuss parabolic equations, but the theory can easily be 
modified to include this case [Obc06]. The focus of [Obe04, Obe05, Obe08a] was on 
numerical approximation. Well posedness for the discrete equations was not estab- 
lished. This result was not needed for convergence, because a small perturbation 
of the equation makes it proper without affecting consistency. 

2. PDES ON WEIGHTED GRAPHS, DEFINITIONS, PROPERTIES AND EXAMPLES 

We consider a finite weighted directed graph G = (V,E,w). Here V is the set 
of vertices {x\, . . . xn} and E C V x V is the set of oriented edges. We denote by 
e xy the edge {x, y) that goes from the vertex x to the vertex y. The positive weight 
function w: E — > R + is defined on directed edges that join different vertices. We 
write 

w(e xy ) = w xy > 

for the value of the weight function on the edge e xy . 

We also identify a non-empty subset dV C V, which we call the boundary of the 
graph. The interior of the graph is the set of vertices V \ dV . We write x, y for 
typical vertices in V. The vertex y ^ x is a neighbor of the vertex x, if there is an 
edge e xy from x to y. The degree of a vertex d{x) is the number of neighbors of the 
vertex x. We shall always assume that for an interior vertex x € V \ dV we have 
d[x) > 1. For each x <E V we fix an ordering of the neighbors of x, which we write 
as 

y(x) = (yi,y2, ■ ■ -,Vd(x))- 

This is merely for notational convenience because our results will be independent 
of the choice of ordering. We also define the set formed by the vertex x and its 
neighbors N(x) = {x,y 1 ,y 2 , . ■ . ,y d ( x )}- 

The directed distance between x and its neighbor y x is 

d(x,y) = —. 
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Figure 1. Regular convex geometric graphs: the boundary con- 
sists of the extremal vertices. 



We also set d(x, x) — 0. The distance between two arbitrary vertices x and y is the 
minimal path distance 



where the infimum is taken over all finite paths proceeding via neighboring vertices 
which start at x and end at y. If there is no path connecting x and y we set 
d(x,y) = +oo. We say x is connected to y if d(x,y) is finite. We say that the 
graph is connected to the boundary, if for any x <E V there is some y £ dV with 
d(x, y) < oo. 

Example 1. We define geometric graphs, such as those pictured in Figure 1. In a 
geometric graph the edge relations are determined by a set of I non-zero linearly 
independent vectors {vi,V2, ...,«;} in K n , where I < n, as follows. Set 



A boundary vertex is a vertex that does not have full set of neighbors. A compact 
geometric graph (as a metric space) is necessarily finite. A convex geometric graph 
is defined by the property that for any x € V, if we have that y = x + mvj 6 V for 
a natural number m, we also have 



A regular convex geometric graph is one that has at least one interior vertex; in 
other words, for some x € V we have x + H C V. The interior of V is the set of 
interior vertices, and the boundary dV is its complement, which consist of precisely 
those vertices which do not have a full set of neighbors. The symmetric weights for 
adjacent vertices are given by inverse of the Euclidean distances between vertices 



We generally use boldface type to emphasize that a quantity is a vector. 
Definition 1 (Vector Order Relations). Given vectors p, q € R" we write 



d(x, y) = inf \ d(x, xi) + d(x 1 ,x 2 ) + ... + d(x k ,y) 



H = (ui,-«i,« 2 , -i>2, -vi)- 
For an interior vertex x, the neighbors of x are given by 

y{x) =x + H. 



x + nvj € V, n = 1, 2, . . . m — 1. 



w(x, X ± Vj) = l/||l>j | 
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Also we write 

p < q 4=>- p < q and pi < q L for some i. 

We also write 1 = (1,1,. ..,1), where the number of entries will be clear from 
context. 

2.1. PDEs on graphs. The class of all functions u: V — > K will be denoted by 
C(G). The tangent space TG(x) to the graph G at an interior vertex x is 

TG(x) = R d( - X \ 

where d(x) is the degree of x. The tangent bundle TG is the disjoint union of the 
tangent spaces 

TG= |J TG{x). 

x£V\dV 

Given a function u € C(G) we write 

u(y( x )) = (u(yi), u(y d{x) )), 

for the list of values of u at the list of neighbors y(a;) at interior vertices x £ V \dV. 
Definition 2 (Gradient). The gradient vector operator 

V: C(G) i y TG 

acting on a function u E C(G) associates an element of TG(x) to each interior 
vertex x as follows 

(2.1) Vu(x) = (w(e xyi ) (u(yi) - u{x)) w{e xyd ){u{y d ) ~ u{x))) . 

We consider a special class of operators 

F : C(G) -> C(G), 

called Partial Differential Equation (PDE) operators. These operators are local in 
the sense that the value of F(u)(x) at an interior vertex x depends only on the 
value of the function u and its gradient Vit at x. 

Definition 3. We say that F : C(G) — » C(G) is a partial differential equation on 
the graph G, if it can be written in the form 



(PDE) 



F(u)(x) = f{x,u{x),Vu(x)) for xE V\dV, 
F(u)(x) = g(x) — u{x) for x € dV, 



where / : (V \ dV) x K x TG -> K and g : dV -> E. 

The Dirichlet problem for our PDE requires finding a function u G C(G) such 
that 



(D) 



/(i,«(i),Vu(e))=0, for a; ey\cH^, 
g{x) — u{x) = 0, for x G 9V. 



Note that for u to be a solution to the Dirichlet problem (D) F(u)(x) = must 
hold for all x € V. 
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3. Examples of PDEs on graphs 

3.1. Specific Examples. Consider a finite weighted directed graph G = (V, E, w). 
In the following examples, where the condition on dV is not explicitly mentioned, 
we assume that x € V \ dV is an interior vertex. 

The most studied PDE on a graph is the graph Laplacian. 

Definition 4. The (weighted) graph Laplacian in G is given by 

d(x) 

(3.1) Au(x) = 1 • Vu(x) = 2J w xyt (u (yi) - u(x)). 

i=i 

We define the component-wise maximum and minimum functions on vectors 
x £ M. d by min(x) = min{i£i, . . . , Xd] and max(x) = maxjxi, . . . , Xd}. This leads 
to two natural PDEs related to the graph distance known as the eikonal operators. 

Definition 5. The positive and negative eikonal operators on the graph are 

(3.2) ||Vu(a;)|| + = max(Vit(x)) = max w xyi (u(yi) — u(x)) 

l<i<d(x) 

and 

(3.3) || Vu(x)|| _ = min(Vu(x)) = min w xyi {u(yi) — u(x)) 

l<i<d(x) 

respectively, where we have used definition (2.1). 

Example 2. Consider the homogeneous Dirichlet problem (D) for the positive eikonal 
operator 

J ||Vu(.x)|| + - 1 = 0, x£V\dV, 

[ -u(x) = 0, x e dV. 

The solution is the negative distance function to the target set dV. The positive 
distance function is the solution of the corresponding Dirichlet problem for the 
negative eikonal operator ||Vii(n;)||~ + 1 = 0. Let us verify the last statement. 
Given u(x) equal to the (positive) distance function. Then the minimal component 
of Vu will be a negative one, along the neighbor which contains the shortest path 
to the target set. Because the path from this neighbor will be one vertex shorter, 
the equation holds. We show uniqueness below in subsection 5.2. 

Definition 6. The infinity Laplacian on the graph is given by 

(3.4, Mi) = ! Mt±JMt. 

Definition 7. The 1-Laplacian operator on the graph is given by 
(3.5) Aiu(x) = median(Vu(a;)). 

The median of the set {x±, . . . ,Xk} is found by arranging all the numbers from 
lowest value to highest value and selecting the middle one. If k is even, the median 
is defined to be the mean of the two middle values. 

This operator is studied in [Obe04], where it is shown to correspond to the 
operator in the equation for motion of level sets by mean curvature, provided the 
neighbors are arranged close to uniformly about a circle of small radius centered at 
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3.2. Definitions: Elliptic PDEs and the Comparison Principle. 

Definition 8 (Comparison principle). Given a graph G and a partial differential 
equation F, we say that the comparison principle holds for F in G if we have 

(Comp) F(u)(x) > F(v)(x) for all x € V =$> u(x) < v(x) for all x € V. 

Uniqueness of solutions clearly follows from comparison, since if u and v are 
solutions F(u)(x) — F(v)(x) — for all x € V, so that u(x) < v(x) and u(x) > 
v(x). Thus we conclude u(x) = v{x) for all x 6 V. 

Remark 2. For the Dirichlet problem (D) the inequality u(a:) < for a; G 9V 
follows immediately from > F(v)(x) for x € so that to establish the 

comparison principle we would need to show that u{x) < v(x) for interior vertices 
xeV\dV. 

The following definition refers to properties of the PDE F on interior vertices 
V \ dV. 

Definition 9. The PDE F is elliptic at the vertex x £ V \ dV if we have 

(3.6) r>s, p < q f(x,r,p) < f(x,s,q). 

for all r,seR and all p, q e R d{x) . 

The PDE F is elliptic in G if it is elliptic for all x € V \ dV. 
The PDE F is proper at the vertex i£l/ \ 9V if 

(3.7) r > s, p < q /(x,r,p) < /(a;,s,q). 

or all r,s£l and all p, q E R d< - X \ 

The PDE F is proper in G if it is proper for all x £ V \ dV. 

The PDE F is uniformly elliptic at the vertex x £ V \ dV if we have 

(3.8) r>s, p < q f(x,r,p) < f(x,s,q). 

for all r,seR and all p, q E R d( ~ x \ 

Recall that strict inequality means strict in at least one component. 

The PDE F is uniformly elliptic in G if it is uniformly elliptic for all x € V\dV. 

Example 3. The graph Laplacian is uniformly elliptic while the trivial equation 
F(u)(x) — g(x) — u(x) for all x € V is proper. The eikonal equation and the 
infinity Laplace equation are elliptic, but they are neither proper nor uniformly 
elliptic. 

Observe that if the function u has a local maximum at the vertex x € V \ dV 
we have 

(3.9) Vu(x) < 0. 

If the local maximum is strict, the condition becomes 

(3.10) Vu(x) < 0. 

Remark 3. The condition (3.9) plays the role of the familiar second derivatives 
condition D 2 u(x) < at a local maximum x of u, when u is twicc-diffcrcntiablc. 

The following lemma is a restatement of Definition 9 in terms of local maxima 
of u — v. 
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Lemma 1. Let x <E V \ dV be an interior vertex. The PDE F is elliptic at the 
vertex x if we have 

(3.11) u(x) > v{x) and Vu(x) < Vv(x) => F{u){x) < F(v)(x). 
The PDE F is proper at the vertex x if 

(3.12) u(x) > v(x) and Vu{x) < Vi>(i) => F(u)(x) < F(v)(x). 
The PDE F is uniformly elliptic at the vertex x if 

(3.13) u(x) > v{x) and Vtt(i) < Vv(x) F{u){x) < F{v)(x). 

The hypothesis of (3.11) is precisely that u — v has a non-negative local maximum 
at the vertex x, while in (3.12) the hypothesis is that u — v has a strictly positive 
local maximum, and in (3.13) the hypothesis is that u — v has a non-negative strict 
local maximum. 

Remark 4. These definitions correspond in the continuous case to the definitions 
of proper, elliptic, and uniformly elliptic equations in [CIL92], with the inequalities 
reversed. The reversal of inequalities comes from the fact that some authors, us 
included, use the convention that Am is elliptic, which others use the convention 
that — Au is elliptic. 

Remark 5. Proving the converse of the inequalities in Lemma 1 leads to uniqueness 
for solutions of the PDEs. We will prove that solutions of uniformly elliptic and 
proper PDEs are unique. We will also give other conditions on (degenerate) elliptic 
schemes which lead to uniqueness. 

3.3. Generalized p-harmonious equations. We consider elliptic PDEs that are 
expressed in terms of functions f(x, p) for interior vertices x and that are indepen- 
dent of u. 

Definition 10. The elliptic function f(x, p) is homogeneous if there exists a non- 
negative function wo : V \ dV —> R such that we have 

(3.14) w (x) min(p) < f(x, p) < w (x) max(p), 
for all x £ V \ dV and p e . 

Definition 11. We say F is a (generalized) p -harmonious equation if it is of the 
form 

(3.15) F(u)(x) = w+{x)\\\7u(x)\\ + + w-{x)\\\7u(x)\\- + f(x, Vu(x)), 

where w + and u>_ are strictly positive functions in V \ dV, and the function / is 
elliptic and homogeneous as in (3.14). 

Example 4. Consider the special case of (3.15) given by 

(3.16) F(u)(x) = w Au(x) + wiAiu(x) + w+\\Vu(x)\\ + + w_||Vu(:e)||~, 

where wq, Wi, w+, and w~ are non- negative functions on interior vertices. We 
call a solution of F(u)(x) = with F given by (3.16) p-harmonious with drift by 
analogy with the continuous case in [MPR12]. Functions of this type arise as the 
approximations of p-harmonic functions [Obel2,MPR12] when w + = u>_. 

The classical (not on a graph) normalized p-harmonic operator for 1 < p < oo is 

A"u = -\Vu\ 2 - p A p u = -|Vm| 2 - p div(|Vu| p - 2 Vn), 
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(3.19) F p (u)(x) = -A lU (x) + — ||Vu(x)|| + + — ||Vu(a:)|| 



while for p = oowe set 

(3.17) A^u = IVul^Aooit = |Vu|~ 2 (D 2 uVu, Vzt). 
It is easy to see that 

(3.18) A»u=±A?u+UZu, 

where q is the Holder conjugate of p, 1/p + 1/q = 1. Note that A^u = (l/2)Au 
and that 

Au = Afu + A^u. 

Therefore if we set ujq = 1, u)\ = 1/p, and w+ = w_ = l/(2q) in (3.16) we obtain 
the discrete analogue of the normalized p-Laplace operator 

1 p A 1 u( X ) + j q \\Vu( X )\\+ + ± 

In the latter case an approximating sequence is generated by running tug-of- 
war stochastic games with noise of decreasing step-size. The value of the game 
function satisfies a nonlinear equation, which is directly linked to the existence and 
uniqueness of the solution to the p-Laplacian for 1 < p < oo, as demonstrated in 
[PSSW09], [PS08] and [MPR12]. 

Definition 12. We say F is a positive eikonal equation if it is of the form 

(3.20) F(u)(x) = w + (x)\\Vu(x)\\ + + f(x, Vu(x)) - g(x), 

where w+ and g are strictly positive functions, and the function / is elliptic and 
homogeneous as in (3.14). 

4. Proofs: uniformly elliptic and proper cases 

The uniqueness proofs for different classes of PDEs share the same structure: (i) 
assume a failure of the comparison principle, (ii) identify a set where the comparison 
principal fails maximally, and (iii) use the properties of the equation to propagate 
this set to the boundary and obtain a contradiction. To this end we give the 
following definition. 

Definition 13. Given u, v G C(G) define the set of vertices where maximum of u — v 
is attained 

W = W(u, v) = argmax(u(i) — v(x)) 
and the actual maximum 

M = M(u, v) = max(u(x) — v(x)). 

4.1. Counterexamples to well-posedness. 

Example 5 (Failure of existence). Existence can fail for the Dirichlet problem for 
an eikonal equation with the wrong sign, e.g. |jVw|| + = — 1. A simple example is 
K 3 , the fully connected graph with three vertices V — {A, B, C}, three edges e^s, 
esc, and ecA, and unit weights. Let C be a boundary point, with u(C) = 0. The 
PDE can be written after a simple manipulation as 

u(A) = max{u(S), 0} + 1, u(B) = max{w(/l), 0} + 1, 

which has no solution, since the first equation implies u(A) > u{B), and the second 
u{B) > u(A). 
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Figure 2. The graph and the boundary values. 



In the continuous setting, the PDE for level set motion by mean curvature is 
well-posed in the parabolic case, but uniqueness fails for the Dirichlet problem, 
see [SZ94]. Uniqueness also fails for the Dirichlet problem in the discrete case, as 
the following example shows. 

Example 6 (Failure of uniqueness). Uniqueness can fail for the 1-Laplacian (3.5). 
Consider the planar geometric graph G illustrated in Figure 2, where the weights 
are set equal to one. The twelve vertices in G are 

V=Ui,j),i,je {-2,-1,1,2}: + <3|. 

We set dV = {(i,j) € V: \i\ + \ j\ = 3} (eight vertices). The edges are given by 
segments parallels to the axis in such a way that each for the four vertices in 
V \ dV has four neighbors. We assign boundary values by setting 

g(l, 2) - g(-l, 2) = g(-l, -2) = g(l, -2) = -1 

and 

5 (2, 1) = g(2, -1) = g(-2, 1) = g(-2, -1) = 1. 
Then, the constant functions 

u(i,j) = i, (i,j)eV\dV 

and 

v(i,j) = -i, (i,j)eV\dV 

are both solutions. This is because in both cases, three of the four neighbors have 
the same value. 

Remark 6. Uniqueness fails for the non-homogeneous infinity Laplacian equation 

Aoo-u = g, 

in the case where g(x) changes sign as noted in [PSSW09]. For failure of uniqueness 
in the continuous case see [AS 12]. 

4.2. Proofs. 

Theorem 1. Consider a graph G which is connected to the boundary dV C V. 
Suppose F is uniformly elliptic, then the comparison principle (Comp) holds for 
the Dirichlet problem (D). 
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Proof. Suppose that F(u)(x) > F(v)(x) for all x £ V . It follows that u(x) < v(x) 
for boundary vertices x £ dV. We wish to show that u{x) < v(x) for interior 
vertices x £ V \ dV. Define the set W and the quantity M as in Definition 13. 
Assume that M > 0. By assumption, there is a path from any point in W to some 
point in dV , so by following the path until the first point outside the set is found, 
there exists x £ W with neighbor y such that y £" W. We conclude that a; is a strict 
positive local maximum of u — v, so that Vu(x) < Vji(i). Since F is uniformly 
elliptic at x, we can apply (3.13) to obtain 

F(u)(x) < F(v)(x), 

which contradicts the uniform ellipticity assumption and completes the proof. □ 

Theorem 2. Consider a graph G which is connected to the boundary dV C V. 
Suppose F is proper, then the comparison principle (Comp) holds for the Dirichlet 
problem (D). 

Proof. The proof is similar to that of previous theorem. In this case we take x £ W, 
a positive local maximum of u — v so that Vu(i) < Vi>(cc). Because u(x) > v(x) 
and F is proper we deduce 

F(u)(x) = f(x, u(x), Vu(as) < f(x, v{x), Vv{x)) = F(v)(x), 

in contradiction with our hypothesis. □ 

Remark 7. We thank the anonymous referee for pointing out that indeed what is 
used in the proofs of the above theorems is the condition 

r > s, p < q f[x,r,p) < f(x,s,q) 

for interior vertices x € V \ dV. This conditions is weaker than both uniform 
ellipticity and properness. 

5. Uniqueness for Eikonal and ^-harmonic Equations 

In contrast to the uniformly elliptic case, where (3.13) holds, for the equations 
considered in this section, a strict local maximum of u — v does not imply a strict 
inequality for the values of the equation. Instead, we need to locate a neighboring 
value which is active (see Lemma 2 below) in the equation. 

5.1. A Lemma on Propagation of Maxima. We begin with the definition of 
an active neighbor. 

Definition 14. Given an interior vertex x and a neighbor z of x, we say z is an 
active neighbor of x if 

||Vu(x)|| + = w xz {u(z) - u{x)). 
Lemma 2. Suppose that F is of the form 
(5.1) F(u)(x) = w + {x)\\Vu{x)\\ + + H(u){x) 

where w+(x) > for all x, and the PDE H is elliptic. Suppose that 

F(u)(x) > F(v)(x) 

for all x £ V\dV. Consider the vertex set W and the quantity M as in Definition 13 
and suppose that M > 0. Then we have that if x £ WC\ (V \ dV) and z is an active 
neighbor of x, then z £ W. 
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Proof. Choose a vertex x G W PI (V \ dV). Because u — v has a local maximum at 
x we have Vu(x) < Vv(x). Therefore, by the ellipticity of H (3.11) we have 

H{u){x) < H(v)(x). 

Using our hypothesis on F we conclude 

||Vh(z)||+ >||V^)|| + , 

and because z is an active neighbor, 

w xz (u(z) - u(x)) = ||Vu(a;)|| + > \\Vv(x)\\+ > w xz (v(z) - v(x)), 

where the last inequality holds by the definition of the Eikonal operator. This last 
result implies 

u(z) — v(z) > u(x) — v(x). 

But from the definition of W this last inequality must be an equality. So z is also 
in W. □ 

5.2. Proofs: uniqueness for eikonal equations. 

Theorem 3. Consider the Dirichlet problem (D) for the positive eikonal equa- 
tion (3.20) on a graph G with connected nonempty boundary dV . Suppose u,v are 
solutions of the Dirichlet problem for F . Then 

u(x) < v(x) for x € dV u{x) < v(x) for x G V. 

Remark 8. A similar result holds for negative eikonal equations. 

Proof. Define M and W as in Definition 13. Suppose M > 0. We also define 

C = maxu(i), Z — {x £ W | u(x) = C}. 

Consider x G Z and choose z to be an active neighbor of x. From Lemma 2, 
since (3.20) can be written in the form (5.1), we conclude that z € W. 
Because it is a solution, we have 

= w + (x)\\Vu(x)\\ + + f(x,Vu) - g(x) 

< w + (x)\\\7u(x)\\ + + w \\\7u(x)\\ + - g(x) 

by property (3.14) of /. Because g is positive we must have 

||Vw(a;)|| + > 0. 

This means u(z) > u(x) = C which is a contradiction to z € W. □ 

5.3. Proof of uniqueness for p-harmonious equations. 

Lemma 3 (Of Harnack type). Let F be p-harmonious (3.15) and x G V \ dV be an 
interior vertex. If F(u)(x) = 0, then either 

(5.2) ||Vu(x)||- < < ||Vu(x)|| + 
or 

Vu(x) = 0. 

Note that (5.2) is equivalent to 

(5.3) min?i(y(x)) < u{x) < maxit(y(a;)), 
which is the reason for describing the result as of Harnack type. 
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Proof. Assume V'u(x) ^ 0. Immediately following from the definition, we conclude 
that one of ||Vu(x)|| _ or ||Vu(x)|| + is nonzero. Suppose that ||Vu(x)||~ 7^ (the 
argument will follow in a similar way for the other case). 
Applying (3.14) to (3.15) we obtain 

(i) (wq(x) + w_(x))||Vit(x)||~ + w + (x)||Vu(x)|| + < 

(ii) < w-(x)\\Vu(x)\\~ + (w (x) +w + (x))\\Vu{x)\\ + . 

If ||Vu(x)||~ > 0, then by dehnition ||Vit(x)|| + > which contradicts (i). So 
that we must have ||Vit(x)|| _ < 0. But then using (ii), we obtain ||Vw(x)|| + > 0, 
thereby obtaining the result (5.2). □ 

Theorem 4. Consider the Dirichlet problem (D) for the p -harmonious function F 
on a graph G which is connected to its nonempty boundary dV. Suppose u,v are 
solutions of the Dirichlet problem for F. Then 

u(x) < v(x) for x € dV =>■ u{x) < v(x) for x € V. 

Proof. Define the vertex set W and the quantity M as in Definition 13. Suppose 
M > 0. Also define 

C = maxu(i), Z — {x £ W \ u(x) = C}. 

Consider x G Z and choose z to be an active neighbor of x. Since the p- 
harmonious equation is also of the form (5.1), we can apply Lemma 2, to conclude 
that zeW. 

Next, we claim that Vu(x) = 0. Since z 6 W, by the definition of C we have 
that u(z) < C. On the other hand, u(x) = C. Since we assumed that u is a 
solution, we can apply Lemma 3, to conclude Vu(x) = 0. This implies that the 
set N(x) C Z C W. It follows that we can find a path from x to dV which stays 
in Z . But we assumed that u < v on dV and that u > v on Z, and hence a 
contradiction. □ 



6. Existence results 

The existence theory is somewhat delicate, since there are known examples where 
existence can fail. In the special case of operators of the form 

-Xu + F{u) 

with F elliptic, existence and uniqueness results were established in Theorem 7 of 
[Obe06]. Therefore in this section we consider equations involving no such depen- 
dence on u. 



6.1. Existence for homogeneous equations. We consider homogeneous equa- 
tions, which include the Laplacian, the Infinity Laplacian, the Eikonal equations, 
and the other examples from subsection 3.1. 

Theorem 5. Let /(x, p) be elliptic (3.6) and homogeneous (3.14). Suppose that f 
is continuous in p for each x. Then the Dirichlet problem (D) has a solution. 

Proof. We use the Brouwer fixed point theorem: a continuous function from a 
convex, compact subset K of Euclidean space to itself has a fixed point. 
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We identify the set of functions on the graphs C(G) with R N and consider the 

set 

K = {ue C{G) | u(x) = g[x),x G dV, m < u(x) <M,xeV\ dV} , 

where 

m = min g(x), M — max g(x). 

x£dV x£dV 

Then if is a convex and compact set. 
Define the constant 

L = max wo(x)w xv , 

x£V\dV,y£y(x) 

and supposed that L > 0. Define the mapping T : C(G) — > C(G) by 

(6.1) T{u){x)= u{x) + ^f(x,Vu{x)), lotx£V\dV 

T(u){x) = g(x), for x € dV. 

We claim T is continuous and that it takes K to K. The result will follow from 
this claim and the Brouwer fixed point theorem. 

The continuity of T follows from the continuity of f(x, p) in the variable p for 
each fixed x G V \ dV. To see this, enumerate the vertices 

{xi,x 2 , ■ ■ ■ ,x n ,x n+1 , . . .x N }, 

where the first n vertices are interior vertices and the last N — n are boundary 
vertices. The map T: M. N i— > K w is given by an expression of the form 

(Tu)i =Ui + ^rf(xi, {Uj - Ui}j), 

where the second argument of / is just Vu(x^) for 1 < i < n and by 

(Tu)i = g(xi), 

which is a constant value in it, for n < i < N . 

Note that from (3.14) we have for interior vertices x 

w (x)\\p\\~ < f(x,p) < w Q (x)\\p\\ + , 

so by the definition of L and the positive eikonal operator (3.2) 

~/(x,Vu(z)) < ^w (x)\\\7u{x)\\ + 

= y!D (i) max w xyi (u(yi) - u(x)) 

L/ l<i<d(x) 

= y w o( x ) ( u (yio) - u(x)), for some i Q e [1, d{x)] 
-w (x) w X y H) (u(y io ) - u(x)) 



< 
< 



, i . ,■ • i i it.. . \ 1 1 \ i i ; • — 1 1 l i i 1 

L 

T w (x) w xyt max (u(y) - u(x)) 



< max {u(y)-u(x)}, 

yeN{x) 

where we have used the set N(x) in the last two lines in case ||Vw(a;)|| + < 0. 
Similarly by using (3.3), we get 

T f(x,Wu(x))> min {u(y)-u(x)}. 
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We conclude that 



min u(y) < T(u)(x) < max u(y), 

yGN(x) yeN(x) 



and so 

m < T(u)(x) < M, for all x £ V. 

Then the claim follows when L > 0. When L = we must have wo(a;) = 0, and 
thus f(x, p) = 0, for all interior vertices x £ V \ dV. In this case, any extension of 
g from dV to V is a solution to the Dirichlet problem (D). □ 

6.2. Existence for eikonal equations. We generalize the fact that the distance 
function to the set dV is the solution of the eikonal equation, Example 2, to con- 
struct solutions of more general eikonal equations. 



Theorem 6. The Dirichlet problem 
f(x,u,Vu(x)) = 



||Vu(x)||+ -h(x), x£V\dV, 
u(x) — g(x), x £ dV. 



for h(x) > 0, has a solution. 



Proof. First notice that when g = and h = 1 for all x, the solution is given by 
the (negative) minimal distance along the graph to the boundary vertices 

u(x) — — min dist(x,y). 

y£dV 

It follows that ||Vw(a;)|| + = 1, since w xy {u(y) — u(x)) — 1 when y, the active 
neighbour in the minimum, is the next point on the minimizing path which achieves 
the distance to the boundary. 

Next, we can assume that g > 0, by replacing u by u + c for an appropriate 
constant. We can eliminate the nonzero Dirichlet values as follows. Suppose x £ 
dV, with g(x) = g\. Then add a single neighbor y\ ^ x so that y(yi) = {x} with 



XIJL 



l/gi- Consider the new equation on the extended graph, where x is now 
an interior vertex and y\ is a boundary vertex, and we set g(yi) = 0. In this way 
we have constructed a new Dirichlet problem for the eikonal equation with zero 
Dirichlet values, which yields a solution of the original problem. 

Finally, for nonconstant h, we can consider the graph with w xy replaced by 
w xy h(x). Although this graph may have non-symmetric weights, but this poses 
no additional difficulties. Now we can apply the previous existence result to the 
rescaled equation. □ 



7. Connections with Finite difference approximations 

This section focuses on how PDEs on graphs can be obtained as elliptic finite 
difference approximations to elliptic partial differential equations. The natural class 
of difference schemes for fully nonlinear elliptic equations is the class of monotone 
(elliptic) finite difference schemes, because they respect the comparison principle. 
Such schemes can be shown to converge to the unique viscosity solution of the 
underlying PDE [BS91]. A systematic method for building monotone difference 
schemes was developed in [Obe06], where the class of elliptic difference schemes 
(which agrees with elliptic PDEs, below) was identified. 
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We build elliptic finite difference approximations of the PDE operators by using 
Taylor approximations. For example, for a smooth function u(x) we have 

u(x + h) = u(x) + hu'(x) + ytt"(a:) + % u "'^ + °( h ^ 

and 

u(x -h) = u{x) - hu'(x) + - Ji u "'( x ) + °( h4 y 

Averaging these equations gives the familiar finite difference expression 

u(x + h) - 2u(x) + u(x - h) 2 
u (x) = — 2 + 0(h ). 

First order expressions for u x are also easily obtained. Combining these gives 
\u x \ = — ma,x(u(x + h) — u(x), u(x — h) — u(x)) + O(h) 



and 



-\ u x\ = — mm(u(x + h) — u(x), u(x — h) — u(x)) + O(h), 
h 



which are monotone schemes for the operators |u x |,— \ux\> respectively [Obe06]. 
Likewise, the approximation for the infinity Laplace operator (3.17) is given by 

(7.1) - A 0o n(y) = \ [2u(x) - min u(x) - max u(x) ) + 0(r), 

which comes from the fact that 

Vit 

ir = arg min u(x) = x — r- + Oir 2 ) 

xedB r (y) |Vu| 

and likewise for the maximum 

y + = arg max u{x) — x + r . + Oir 2 ). 
° x£dB r ( y ) |Vu| 

Thus, we can write 

r 2 

u(y+) = u(x) + r\Vu(x)\ + — , 9 (D 2 u(x)\7 u(x) , Vu(x)) + 0(r 3 ) 

I Vur 



and likewise 



r 2 

u(y~) = u(x) - r\Vu(x)\ + -. ^(D 2 u(xWu(x), Vu(x)) + 0(r 3 ) 

so that we have 

«&±«£} _ u[x) = \ (DMx)^, WL) r 2 + (r») 
2 2 \Vu(x)\ |Vu(a;)| 

which leads to the expression (3.17). 

The discretization of the expression onto a grid is less accurate, because of the 
lack of directional resolution, but using a wide stencil grid, a consistent, convergent 
scheme can be obtained [Obe05]. In fact, the expression (7.1) is second order 
accurate for smooth functions [Obel2]. 
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Another example is the equation for the convex envelope [Obe07]. The finite 
difference schemes for this equation lead to a PDE on a graph [Obc08b, Obe08a]. 
Neither the Obstacle problem for the convex envelope, nor the Dirichlet problem for 
the convex envelope [OS11] are proper. The PDE has the form F(D 2 u) — Xi(D 2 u), 
where Ai is the smallest eigenvalue of the Hessian D 2 u. This is discretized using 

, , . u(x + hv) — 2u(x) + u(x — hv) 
X 1 (u)(x)= mm '- -V -+0{h 2 ). 

\\v\\=i hr 

Consistency follows from the fact that the smallest eigenvalue of a symmetric matrix 
M is given by min||„|| =1 (Af • v,v). Applying this fact to the Hessian matrix, and 
using the Taylor series computation gives the result. 

8. Conclusions 

We have proved uniqueness and existence results for a wide class of uniformly 
elliptic and degenerate elliptic PDEs on graphs. The approach combined viscosity 
solutions type techniques (the comparison principle) with connectivity properties of 
the graph to establish uniqueness results. Existence results were established using 
fixed point theorems, or, in the case of generalized distance functions, explicit 
solutions. 

In the future, we hope to solve the discrete equations on unstructured graphs. 
In particular, we hope to implement fast solvers. 

In addition the existence theory could be extended: natural analogues of the 
existence theory for viscosity solutions (barriers, Perron's method) could be studied. 
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